Gamma Rays from Heavy Neutralino Dark Matter 
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We consider the gamma-ray spectrum from neutralino dark matter anniliilations and sliow that 
internal bremsstrahlung of W pair final states gives a previously neglected source of photons at 
energies near the mass of the neutralino. For masses larger than about 1 TeV, and for present 
day detector resolutions, this results in a characteristic signal that may dominate not only over the 
continuous spectrum from W fragmentation, but also over the 77 and 7Z line signals which are 
known to give large rates for heavy neutralinos. Observational prospects thus seem promising. 
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INTRODUCTION 

One of the most favoured dark matter candidates since 
20 years or so, is the lightest supersymmetric parti- 
cle 0, 0. Most of the expected phenomenology has 
been worked out, and is implemented in freely avail- 
able, extensive computer packages like DarkSusy 3] 
and MICROMEGAs . The vast majority of the analy- 
ses has been performed in various constrained versions of 
the minimal supersymmetric extension of the Standard 
Model (MSSM) (for reviews, see e.g. Q), where either 
radiative breaking of supersymmetry or the GUT condi- 
tion on gaugino masses (or both) have been imposed. It is 
clear from all these studies that a suitable neutralino can- 
didate for dark matter is indeed available, with relic den- 
sity as measured by WMAP (f^cDM^i^ = 0.113 ± 0.009, 
where ricDM denotes the ratio of cold dark matter to 
critical density and h is the Hubble constant in units of 
100 kms"^ Mpc^^) 6], and with great prospects of being 
detected either at the CERN LHC, or at the various di- 
rect and indirect detection experiments of halo dark mat- 
ter that are currently operational or being constructed. 

However, one may also ask the unpleasant but not un- 
realistic question what happens if supersymmetry is not 
found at the LHC. If supersymmetry is still present, that 
would probably mean that the mass scale of the light- 
est supersymmetric particles is beyond the kinematical 
reach of the accelerator. For neutralino masses at the 
TeV scale, also the scattering cross section for direct de- 
tection would necessarily be small, as would many of the 
rates (antiprotons, positrons) for indirect detection. An 
exception seems to be gamma-ray detection, firstly be- 
cause rates do not fall off as rapidly 7] and secondly 
due to eminent new gamma-ray telescopes, most clearly 
demonstrated by the recent spectacular performance of 
HESS, in particular as regards the multi-TeV signal that 
has been observed towards the galactic center y. With 
such new astrophysical gamma-ray instruments of un- 
precedented size and energy resolution either in operation 
Q ^ or under construction 0, ^| , it is appropri- 
ate to investigate possible levels of signals and spectral 
signatures of heavy dark matter particle annihilation. 



Most of the previous calculations have been carried 
out at tree-level, with radiative corrections typically (and 
correctly) believed to be at the few percent level. In some 
cases, radiative corrections may on the contrary relieve 
the annihilation rate from inhibiting factors having to do 
with the Majorana fermion property of neutralinos and 
the fact that annihilation in the dark matter halo effec- 
tively takes place at rest 1]. One example of this is the 
radiative "correction" of XX ~^ ff through the emission 
of a photon in the final state. Here, the first-order cor- 
rected cross section can be many orders of magnitude 
larger than the tree- level result |l3l |. 

A second example of an unexpectedly large cross sec- 
tion is that of the second order, loop-induced 77 and Z-j 
annihilation 0, where in the high mass, pure 

higgsino (or wino |l6|) limit the branching ratio normal- 
ized to the lowest order rate can reach percent level, 
despite the naive expectation of being 2 to 3 orders of 
magnitude smaller. The origin of this enhancement has 
only recently been fully understood, and is explained by 
nonperturbative, binding energy effects in the special sit- 
uation of having very small (i.e. galactic) velocities and 
very large dark matter masses as well as small mass dif- 
ferences between the neutralino and the lightest chargino 

In this Letter, we focus on gamma rays from neutralino 
annihilation into charged gauge boson pairs and show 
that there is yet another, previously neglected enhance- 
ment mechanism, appearing already at first order in acm- 
radiative processes containing one photon in addition to 
the weak bosons in the final state, give a new source 
of photons which peaks near the highest possible energy 
(the mass of the neutralino). This turns out to be a very 
beneficial effect for the potential detection of neutralinos 
with masses m-^ > 1 TeV, as the normally soft spectrum 
of continuous photons comin g fr om the fragmentation of 
W or Z bosons (see, e.g., [tI 111! FlSj) gets a high-energy 
supplement with a clearly distinguishable signature. This 
is reminiscent of the case of Kaluza-Klein dark matter, 
where internal bremsstrahlung in annihilation processes 
with charged lepton final states dominates the gamma- 
ray spectrum at the highest energies (see also H^l). 
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GAMMA RAYS FROM NEUTRALINO 
ANNIHILATIONS 



In most models, the lightest stable supersymmetric 
particle is the lightest neutralino, henceforth just "the 
neutralino" , which is a linear combination of the super- 
partners of the gauge and Higgs fields, 

X = X? - N^iB + Ni2W^ + 7Vi3#? + NuH° . (1) 

In order not to overdose the universe, a TeV-scale neu- 
tralino must generally have a very large higgsino frac- 
tion, Zh = \N13f + \N14f , if the usual GUT condition 
Ml ~ M2/2 is imposed; otherwise a heavy wino would 
also be acceptable. In the following, we therefore focus on 
higg sino- like neutralinos, with Zh ^ 1 and N13 « ±Nii 
[23. For the high masses we are interested in, the anni- 
hilation rate into charged gauge bosons often dominates. 
Internal bremsstrahlung in these final states are there- 
fore of great interest to investigate. Moreover, we note 
that all our results are almost independent of the relative 
velocity v of the annihilating neutralino pair. Analytical 
expressions are therefore presented in the limit of vanish- 
ing velocity, but should be applicable both at the time of 
freeze-out (w/c ~ 1/6) and to annihilating neutralinos in 
the galactic halo today {v/c ~ 10^'^). 
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FIG. 1: Contributions to XX W'^W 7 for a pure higgsino- 
like neutralino (crossing fermion lines are not shown). 

For a pure higgsino, the only contribution to the low- 
est order annihilation cross section into charged gauge 
bosons comes from a t-channel exchange of a chargino; it 
is given by 
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where m ± and mw are the lightest chargino and W 
masses, respectively. 

Let us now consider radiative corrections with a photon 
in the final state in addition to the W pair. Just as at 
lowest order, the potential s-channel exchanges of Z and 
Higgs bosons vanish, and the only Feynman diagrams 
that contribute are shown in Fig. ^ To zeroth order in 
e = mw/m^, and retaining a leading logarithmic term, 
the resulting photon multiplicity is given by 
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where x = E-Juriy. and 5 = (na^± — m^)/mw- 

Several interesting features can be identified in this 
expression. For large mass shifts 5 the last two terms 
dominate. They originate from longitudinally polarized 
charged gauge bosons in the final state, which are forbid- 
den in the lowest order process because of the different 
CP properties of the initial and final state [Ul . Remem- 
ber that in the limit of vanishing relative velocity, the ini- 
tial state must be an S'-wave with pseudoscalar quantum 
numbers due to the Majorana nature of the neutralino. 
The emission of a photon, on the other hand, will open up 
this channel in the ^Sq partial wave, potentially leading 
to very large cross sections |27j . However, in supersym- 
metric scenarios with a heavy higgsino- like neutralino one 
usually expects a mass shift 5 < e l2ll | , in which case 
the longitudinal part (last two terms) can be neglected. 

For small mass shifts, the cross section is instead dom- 
inated by the production of transversely polarized gauge 
bosons. This results in a peak in the spectrum at high 
energies that becomes more and more pronounced for 
higher neutralino masses. The appearance of this peak 
can be understood by observing that for very heavy neu- 
tralino masses the transversely polarized W bosons can 
be treated as light and thus behave in the same way as 
infrared photons radiated from the neutralino/chargino 
line in Fig. The mechanism that takes place is, in other 
words, an amusing reflection of QED infrared behaviour 
also for W bosons: The kinematical situation when the 
photon and one of the Ws leave the annihilation point 
each with maximal energy, gets an enhancement, since it 
is automatically accompanied by a very soft W . This is 
also reflected in the symmetric appearance of the x 
and a; ^ 1 poles in the first terms of Eq. JSJ. 

As an illustrative example, we have chosen a typi- 
cal higgsino-like MSSM model, fulfilling all experimental 
constraints, as specified in Table (similar models are 
found in, e.g., the focus point region of mSUGRA). The 
resulting photon spectrum from internal bremsstrahlung 
of W pair final states is shown in Fig. 13 The symme- 
try around x ~ 0.5 in the spectrum indicates the related 
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FIG. 2: The photon multiphcity for the radiative processes 
XX W'^W'^. The dots represent the MSSM model of 
Table as computed with the FormCalc package |23l for a 
relative neutralino velocity of 10~^. The thick solid line shows 
the full analytical result for the pure higgsino limit of the 
same model but with zero relative neutralino velocity. The 
thin solid line is the corresponding approximation as given in 
Eq.Q- Also shown, as dashed and dotted lines, are two pure 
higgsino models with a lightest neutralino (chargino) mass of 
10 TeV (10 TeV) and 1.5 TeV (2.5 TeV), respectively. 
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TABLE I: MSSM parameters for the example model shown 
in Fig. 1214) and the resulting neutralino mass (m^), chargino 

mass (m ±), higgsino fraction (Zh), branching ratio into W 

^1 

pairs (W^) and neutralino relic density (Q^h^), as calculated 
with DarkSusy j^] and micrOMEGAs Masses are given 
in units of TeV. 

nature of the peak and the infrared divergence. For com- 
pleteness, we have also included a very high mass (10 
TeV) hi ggsi no model which has received some attention 
recently |l7ll23| (even though thermal production of such 
a neutralino in general gives a too large rJcDMj unless one 
allows for finetuning of parameters like the psedudoscalar 
Higgs mass H^D- In addition, the case of a hypothetical 
model with a very large mass shift is shown (where the 
contributions from longitudinal W bosons dominate at 
high energies). 

Let us now consider those contributions to the gamma- 
ray spectrum from the decay of heavy neutralinos that 
have been studied earlier. Secondary gamma rays are 
produced in the fragmentation of the W pairs, mainly 
through the decay of neutral pions. In addition to the 
secondary spectrum, there are line signals from the direct 
annihilation of a neutralino pair into 77 and Z'y |lfJ | . 
Due to the high mass of the neutralino, these lines cannot 
be resolved but effectively add to each other at an energy 
equal to the neutralino mass. 
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FIG. 3: The total differential photon distribution from xx 
annihilations (soUd Une) for the MSSM model of Table|T] Also 
shown separately is the contribution from radiative processes 
XX ~* W'^W~"f (dashed), and the W fragmentation together 
with the XX ~^ 77) -^7 lines (dotted). 
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FIG. 4: The same spectra as in Fig. |3] as seen by a detector 
with an energy resolution of 15 percent. 

For comparison, again using the model of Table HI 
Fig . 13 shows the contributions from secondary photons 
[l7| and the line signals, as well as the new source of 
photons from the internal brcmsstrahlung diagrams of 
Fig. [2 

The practical importance of the latter contribution can 
be appreciated even more, when considering a finite de- 
tector resolution of 15 %, which is typical for atmospheric 
Cherenkov telescopes in that energy range; the result is a 
smeared spectrum as shown in Fig. 0] One can see that, 
although the strength of the 77 and Z'y lines already are 
surprisingly large [3, the contribution from the internal 
bremsstrahlung further enhances this peak by a factor of 
2. The signal is also dramatically increased at lower en- 
ergies, thereby filling out the "dip" just below the peak; 
this latter effect will of course become even more pro- 
nounced for better detector resolutions. 
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CONCLUSIONS AND DISCUSSION 

In this Letter, we have presented important radiative 
corrections to the gamma-ray spectrum from heavy neu- 
trahno annihilations. They contribute a characteristic 
peak shape at the highest energies, competing with the 
77 and Z7 hue signals in today's detectors. 

When it comes to predicting absolute gamma-ray 
fluxes, there is, unfortunately, a great uncertainty which 
is primarily due to the unknown dark matter cluster- 
ing properties of the Milky Way halo. Particle physics 
alone predicts, in lowest order perturbation theory, a 
flux which falls roughly as m^^. On the other hand, 
the possible nonperturbative enhancement discussed in 
|l7| is important for heavy neutralinos, and may give a 
substantial boost to the gamma-ray signal. The process 
discussed here should benefit from a similar boost as the 
line signal treated in [T^ . since the enhancement is due 
to binding effects in the initial state. Thus, even if the 
radiated W boson is soft compared to the initial state 
total mass, it still causes a virtuality of one of the neu- 
tralinos which is much greater than the binding energy. 
Therefore, the factorization of the process into long dis- 
tance and short distance kernels performed in 17\ should 
be valid also here, and our curves for the ratios of cross 
sections (dA^^ / dx) should not be affected much. 

Since the absolute gamma-ray flux - although not un- 
likely considerable in size - is hard to predict, it is rather 
the spectral shape that eventually may separate a dark 
matter signal from the background. The HESS observa- 
tions of the galactic center, e.g., show a power-law en- 
ergy spectrum with a spectral index of about —2 up to 
at least 10 TeV 0, and for Higgsinos as heavy as that, 
we expect a spectrum that is too hard to explain the 
full data set |2J|. Nevertheless, even a lighter higgsino 
could still partly contribute; once one has access to bet- 
ter statistics, a characteristic distortion in the spectrum 
would then be distinguishable at the dark matter parti- 
cle's mass. One should furthermore bear in mind that 
the best prospects for detection might therefore not be 
found near the galactic center, but rather for sources with 
a low or at least well understood background. Examples 
for this could be nearby dark matter clumps or interme- 
diate mass black holes [l^; a thorough analysis of the 
detectional prospects for these candidates, however, is 
beyond the scope of this Letter. 

Finally, we note that the radiative corrections pre- 
sented in this Letter add to the total annihilation cross 
section, quite independently of the relative velocity of the 
annihilating neutralinos. This is relevant for any precise 
calculation of neutralino relic densities and is therefore 
of relevance for computer packages like DarkSusy and 
micrOMEGAs. 
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